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ABSTRACT 


Mme s or metncmcOnerol! laws Of a bank-to=turn mis- 
Sile using and optimal estimator in the terminal guidance 
phase Were designed, and the effect of increasing the number 
of measurement sensors in the missile to generate more 
mrormatilon On the State was investigated. In the design of 
mie coOncrol law the modern optimal control theory with the a 
quadratic performance index was used. Implementation of 
mes control taw required the use of a Kalman filter as the 
Seeimal G€stimator. Ihe extended Kalman filter algorithm was 
utilized in the present study since the measurement states 
ZaaemiOn-lintcam rumections Of the state vectors. In order to 
Mewemiiic Crttects Of the implementation of the increased mea- 
Surement sensors, two-, four-, and six-measurement sensors 
meme assumed to be implemented in the optimal estimator. By 
computer analysis, the designed guidance laws were evaluated 
and the effect of the implementation of increased measurement 
BemsOrs was tested. 

ieee omor = ric Simtracion revealed that the designed 
memaance law was Successful within the specified scenarios, 
the effect of the implementation of increased measurement 
a iseeOrs fOr the estimator was favorable only in that 
mierecased measurement sensors generated more information 


meout the state vectors, 
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Pe JeNT RODUCT TON 


Miemeane-ce-turn mMisstle has high lift acceleration in 
a direction perpendicular to its wings. For airbreathing 
Missiles which are required for large stand-off ranges it 
Meso eOLters the advantage of lower inlet angles of attack 
Semen that of skid-to-turn missiles. 

Although the control laws using the modern guidance 
Memeno! theory are more complex than the normally used 
Proportional navigation air-to-ground method, they have 
Meet) potential for maneuvering targets in air-to-air 
engagement situations. The application of optimal control 
[meme stimation theory to bank-to-turn missile configurations 
Meee ecieUritca im ecrdér to obtain increased performance. In 
miemcapplication of the optimal control theory to the bank- 
Meeeetirn Missile, it 1S necessary to have information on all 
Peeates Or an Estimate of the state variables that have not 
Meeu measured. Since the state information available from 
miemuypical missile sensors is limited, it is necessary to 
Salpeoy an €Stimator. the estimator used in this study was 
mers t Order extended Kalman filter. 

Mine present work addressed the design and evaluation of 
fmeewoptimal state estimator and optimal control laws for 
Pier cation to a bank-to-turn missile. In the development 


Seeeemis work, Chapter 2 will describe the kinematics of the 


5 





missile and formulate the equations of the motion. All 
for assumptions adreslisted. tJhe optimal controller is 
Peeeenoce ii) Giapter soe) the Simulation results of the 
Raaenolber Ommtaree Sceeiari1os described later are used to 
mie wene sSuUltability of the optimal controller. fhe 
Optimal estimator is covered in Chapter 4, the formulation 
Beechne State equation, Measurement equations and the 

Pee aerOls Ot Cem clemenuts of Jacobean matrix are covered. 
ivemextendcea Kalman filter aloorithm for non-linear system 
Memethen reviewed. Before the final computer simulation a 
discussion is given on the nominal parameters needed in 
Mente nomeye halman filter. Ihe results of the simula- 
Mmiemeor E€StimatOrs With two-; feur-, and six-measurement 
Meetors On a typical scenario are then presented. In 
Seino. tie Simulation results of the control laws imple- 
mented the estimator with two-, four- and six-measurement 
mons wOMmecnree SGenarios 1S described. The detail analy- 
re Ot the results is developed in order to investigate the 
mieecet Of the key variables of the control system on a bank- 
Mee turin missile. the mean miss distance determined from 50 
Monti Carlo runs as a performance standard was used in the 
mirmysis Gr the results of the control laws as a function of 
Maewkey Variables. Finally, the conclusion are summarized 


meee ne last Chapter. 
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ie Cte eee ee ba NA-TO-TURN SMISSILE 


Eee AOS UMPTIONS 


piveeccOMeceny sallam pie weqUations of motion of the missile 


Will be developed, representing the positions in space, 


mader the following assumptions: 


ih 


Diem Oe dimoccaragntemeoordinate svstem will be used 
Acme eceelicuelal he rerence . 

Toms ieee, Che eQuattemeof motion the stability 
Me mates deep eccderorwence missile body. 

(icmiwerwerc VeUOGiIy duc LO thrust is constant. 
ise meommomitoothemacceleration in velocity due to 
light St. 

PoOGhmesencrol Stimace er Surfaces is rigid. 

ceviimie Cm uOmencmmOdveaxesS,,cach Control system has 
Cm mOncm@etreouo, treegom. The missile's control 
lee Vctdemolmyeccor adcussnorman to the velocity 
eoCmetiteeiscesetine dot product Of acceleration 
VCGuOUromeum—mcrocity Vector 1s zero, also’ the 

oe ronaewoneaeramenrouch the missile center of 

Cea by. 


tie sitkerwancles 4 6 and @o will be used to 


Veo 


eeSecripec MEmew Orlemmation of the missile with 


Eeovcel COmiMeteETa mes pacce, where >, and o., are 


horizontal and vertical flight path angles, and 6 


Pmencmuoll Or Dank angle. 
1S 





ae ews eee rece aosSewith time constants i. and 
as fem TOmmeiediee Command Ee ama sa Mem mee 
MoeCelertat fom | Cc Onmana a fede, pe Ole Lede rece: OT 


the dynamic response of the missile. 


meee SiISSTLE CONTROL METHOD 

PiepOre eO1ne gimeo the mathematical detail concerning 
Mime mOtiOn of a missile in space as a result of a guidance 
Pend, tteis helpiul £0 review with the control law for a 
miitestO- turn missile, the guidance system detects whether 
Mitemmicsiice 1s flying too much or too high to the left, 
Mme Or Vertical. the guidance system measures these 
Mameteloms | Or Crrors and sends signals to the control svstem 
Mmmeeieomene se weGrOrs to aero. the task of the control 
meee tmeretoOre 1S tO maneuver the missile quickly and 
meerctently as a result of these signals. In a bank-to-turn 
Pet eeleme seems tie Guidance angular error detector produces 
meemotcnals KR and » in terms of polar coordinate expres- 
mem, showed in Figure 2.1. The same signals can be 
Meee oscasim another way, that 1S, in cartesian coordinate. 
The usual method is to regard the $ signal as a command 
memrort chrough an angle of measured from the vertical 
p@eetien tO Maneuver outwards by means of the missile's 
meeevators. the method of maneuvering the missile is as 
momows. Ihe command goes as a positive command to one 


Bemerol Surface and a negative command to the other, this 


eo 









ESSE 


—ne 


eee ol eeetiiool Le (@ONTR@I, SYSTEM IN POLAR CO-ORD. 


Meetoes the missile to roll. The R command goes to 
Mimraces always aS a positive demand, this causes the 
Memeestte tO accelerate normal to the velocity vector. The 
Mermemction 25 €O make the response in roll fast so that the 
commands can be applied simultaneously which makes for 

Peo icity, Only a pair of control surfaces are used as 
eeecoroens and elevators at the same time, control is obtained 


Peeieans Of the separated servos. 


meee GEOMETRY OF MOTION IN SPACE 

Pim mmoim Grea mss tle as a particle may be described 
Meeusing coordinate measured with moving axes (relative- 
fietrwen analysis). the analysis of motion can be simplified 
ieeeusing measurements made with respect to a moving coordi- 


meee system. In present work, the equations of motion will 


If! 





MISS LEE “NORMAL 
PORCE SAALS 





iE 


MSS iLEE 
LONGITUDINAL 
AXIS 


JmeeitoS ILE WING AXIS 


PUuGUREecnc GEOMETRY OF RELATIVE MOTION. 
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be described with respect BOmemewiNitial missile body axes 
meleonmease tie inertial axes. After launching, the missile 
meayeaxes Changes position, but the relative position, 
Mmemocity and acceleration are still computed in inertial 
mies tat 15, im the imitial body axes. The current zero 
mesOnrtemlss distance 16 Calculated. This will then be 
transformed to values in missile body axes to compute the 
mount Of the Control inputS to the missile using Euler's 
transformation. 

eee OWSeene MisSilewand target aS points with 
Mmespective vector velocity Vi and Ve: iM thes analysis, 
because |V,| is assumed to be at least AS , and the 
em@ete Of attack is assumed to be small, the lead angle 49 


meamvoem.) vv, and the missile longitudinal axis is small. 


R 
mame voevector denote the relative position of the target 
Mmemerespect to the missile. The orientation of the sight 
Memgemvector in inertial Space can De represented by the 
angle YR and op Mosel Merle eboure 2.2.2. [he sieght 
me vector 1s resolved into three components in inertial 


peees as follows: 


Re = R*cos (6p) *ccs (yp) = Xt - Xm eae 
Ry = R*cos($,)*sin(vp) =aye Yn (2552) 
Rz = - R*¥sin ($n) = Zt - Zm 2s) 


fier e R is the magnitude of the sight line vector. 


io 





Dertince the relative velocity vector: 


Where Vv, and V, are target and missile velocities. 
in te oOenateemoretie sloht line vector's magnitude 
Mmeecoual to the component of VR along ke. lt Gan) be 


Beepressed as following: 








2) TO BY ae Zz ee 
a + RR f Ko 
Where ee = Vey = is 
Vey Vey ; my 
eee) mz 


Mmeeeec relatrve velocity components in the inertial frame. 
Mimeiiomeconctry, the relative angles vp and Ip also 


feeeeenic expressed by the relative vectors: 
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Where OR and Vp are Gmemetevation and the azimuth angles 
Seetne relative Sight line vector in the inertial frame. 


Then the rate of change of both angles are 








V 
me to 
R R 
roam COS tev = Smee -Sinw ¥ Poe Osie. Vo. 
: Vv ex Vv ee Ve 1 ¢ 
: ? 2 ae Se) 
Geek + Re”) 
x Me e 
‘ ee, Se cosy. V, 
<< <5 > 2 Ce 
Rcos6. (R~ + R eee 
) i et ee) COs Oo 
V 
niere Vpg and Vey Miecmieme sabes OF Ciance in 6 and 
; v 


Y,, components, Pome nenol~e 2,949 1Mese quantities 
feeebe used later in developing the estimator state 
wyeations and wiil also be used to represent the variables 
Memechie missile seeker. The quantities Op ; - : Vp ; Ye : 
Ee and R Ueemrupeca by Morse) are the set of measurements 


fmmed available from the missile seeker. 


Sees tC NEMATICS OF CONTROL VECTOR 
puolineoe2-4 ard 2.5 are views from hear and side of 


missile. a Wado me Gcimcseienuine Magnitudes of the 
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FIGURE 2.5 RELATIVE VELOCITY COMPONENTS IN POLAR CO-ORD. 
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feemle acceleration due to tne acceleration command (ac) 
and the bank angle due to the roll rate command (Pc) each 
meer. & represents the gravitational effect at the center 
meer avity Of the misSile. As mentioned in control method, 
m@emolgn Of a is such that it 1S positive upward with 
meecoct tO the wings. Figure 2.5 shows that the accelera- 
mueenormal tO the velocity vector has the components due 
Memcne ClevatoOr control input and the effect of the gravity. 
mine tOtal magnitude of the acceleration normal to the 


merocity vector can be represented as 


a *cosd - a*sinée 
m ? = V 


Miemeclocity and acceleration components for simple 
Mipeuiar motion can be represented as 


at 
2. ey 0 


eeewnererence 1. 
same tmis definition, the ditferential equations 
Me erioine the rate of change of the flight path angles are 
a jcOs> - gcoss,, 


i ee ee o 
oe V (Zee 0:) 
mM 
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—— eeu s 
V 


Nee am*cos®@ - g* cose mec Mmacmneudce Of the normal 
acceleration component to the velocity vector in missile 
smoe plane and am*sing@ is the magnitude of the normal 
aa oretton COmponent to tne velocity vector in the x-y 
Meee. the acceleration component of the missile velocity 
Yn ewe ecw eolpolc@s Of thenecttect Of gravity in the 
Meesiile side plane under the assumption that the missile 


Mmemocity 1S constant. The differential equation is: 


V, = -g*siné, eae 


Mienrates Of Change of the missile position components 
(Am, Ym, Zm) can be obtained by integrating the components 
Seemissile velocity. From Figure 2.2, the components of 


Meeoriec Velocity (Vmx, Vmy, Vmz) in inertial frame are 


given by 
y = i = ee * ? as 
x, ieee ie cosé,, cosy. (Zest Ss) 
Y =YV Se eos 6. Sanh OE elles) 
m my m V V 
ew ae a Saal e: (eon 
m Mz m \ 


Roeper emrcspolse Of the Missile tO input command in 
meemal acceleration and roll rate be considered. If a 


@emand is made on a missile for a transverse acceleration, 


a5 





memo miiltiated by Sending a signal to the appropriate con- 
@eeieestriace. The missztle acceleration abe follows the 
femanded acceleration ay in a manner which may be 
Gmaracterized by a frequency wa Pvedemer Cock trequency) 
mmeea damping factor P . For simplicity, in present 

fem, Only the first order lags with time constant Hs and 
T Memo suiled ene wcduaeloOn Of the first-order system can 


O 


Merrepresented as follows: 
X(t) + 1./t#X(t)) = £(t) 


ierone the Laplace transformation on has 
Gop ee ennte eGo ler) | ) 


meee tTO interal Condition. Applying these to the BTT 


feeec ile comtrol system. 


Gea ater= a oes) 


Lee oe nay 


Mearranging the differential equation. 


a = CB - aes fez 8) 
Cee yea es 
eo =P Ce ea 





nhere a. and De O@ew tae acceleration and roil rate 
Semmanas. Equation (Z.20) 21S by definition. From the 
Meee Gitrerential equations the control ratios are 


miatistormed as follows: 


Bese 1 
= a ae) 
a(S) a or 
e a 
| (22) 
PAS.) teoe 1 / 
C Pp 
mepactons (2.10) through (2.22) constitutes the dynamic 
pemiations describing missile motion in response to the 
memimand inputs ac and Pc . The solution of these 
SaiaetOns provide miss1le position (Xm , Ym , Zm) , 
Smmentation and magnitude of velocity Cv, : 8. ania 
meteitatiOn and magnitude of control acceleration (¢ , am) 


E. CO-ORDINATE TRANS FORMATION 

CmcMmutare tne Missile 1S On guided flisht to compute 
Meemmacnitudes of the control inputs at any moment, it is 
Memecsaty tO transtorm the instantaneous values from the 
mere rs in the inertial frame to the instantaneous values in 
meay trame. the transformation from one frame to another 
Saeoc accomplished through a transformation matrix. The 
means tOrmation matrix wiill be developed. 

Define an inertial coordinate system with unit vectors 


Meek, Also,define a missile body axis coordinate system 
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Meer envectors by, }, ke tt is desired to find the trans- 
mermacton between the 1, J, K system and the i, 3, k system. 
ee Ceeiimettewcncnought to be oriented so 
Meee pOINts north, J east, and k down. Similarly, the 

fue System has i alone the longitudinal axis of the 
Meee (which by assumption is along the velocity vector 


fee) Out the right wing, and k down. 


Memsider three successive rotations Wes ; a, and > 
eae v,, MCT etMmeicmibenmcmtime tmertial z axis, and 
meemsorOrms to an intermediate axis systems ly> Jy Ky. 


The O.. CocdrI@i 15 about the 34 oowand Enanstorms to an 
Memeo otem «614, J,, kK. the last rotation 9 about the 
jeects transtorms from flight path axes to missile body 


meee AS a result, the total transformation from inertial 


Memeody axes is given by 


i 1 | ; 
; : eile Jey, 2 J = et J Gre) 


meow 6 61Ss:«[U6trhe «€6tOtal transformation matrix from inertial to 


[oeeeitrame. the elements of the a matrix from Reference Il 


are: 
= : ¥ 
all cos(6,) cos(p,) 
ae * 1 
a1 2 cos(6,) sin(yp,) 
als = 


emia) 
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= COstewioin a) 
a cos(¢) "costs, 


P= sini) “Cams - 


= cos(p)*sini3 )*cos(y) + sin(>)*sin(y) 


eS ee Sint) sin(¢)*cos(y,,) 


cos(9)*cos(s,,) 
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(eee eee Liesl, CONTROLLER 


Peewee ChlowyOnmemem primarily concerned with the 
effect of an optimal estimator on the performance of a bank- 
feecurnm missile, we will first outline the theory of an 
meeima! Controller in this section. The guidance laws 
Mememaetred here were first developed by Stallard [Ref. 10]. 


Meee rtOoliowing is @ brief outline of his work. 


pee )6GEOMETRY 
Pisce ie siOvs ene applicable geometry representing 
meme royected-2ero-Control miss (PZC miss) distance in the 
meet state. The problem is considered to be three- 
ermensional and two major axes system are used: 
H. A Seeker-oriented axis system for estimation or 
measurement, and 
wee tcmelmee Wrincipal axés of the missile for the 
COMmere ly orob lem, 
miperb, 4b denote the target coordinates relative to the 
Meeeeerco atone 1S principal axes. The body system is repre- 
memeed a time TF = 0 by a set of Eulian angles which 
meemte tt tO an Carth fixed SsyStem. The angle Ad is the 
ma@eremental roll angle from the present missile axes to any 
Mire Orientation at time T , shown dashed in Figure 5.1. 
Meme ceemMlss tepresencts the total Projected-Zero effort - 


Mreedistance and 1s defined as the miss distance which 





oecsOceur it no further control was exercised on the 
Beosittewat T= Ti. It is composed of a Y-components 
PP2zCyY miss) and a Z-component (PZCZ [iess|). Under the 
assumption of constant velocity the Peeomponent (PlCh miss) 
meecepresented by the distance that the missile has to ial, 
at any time, since it is not effected bemecontrol, it is 


fee shown in Figure 35.1. 


ioe tame. bE 





PlC-miss 
i ak cCime & 





Memeo PO OMPONENTS OF PeC, MISS DISTANCE. 


oa 





Poet oe DESCRIPTION OF AN OPTIMAL CONTROLLER 


eeeemwea time-Varyine linear system of the forn, 
See le Ce 1 ee at) ea 


Mitere — 15 N-component state vector 
Domi cOmuomentacontrol vector: 
An optimal control is one that minimizes a performance 
mimex Made up Of a quadratic form in the terminal state 
feiss an integral form of quadratic forms in the control and 


ime State: 


de 
1 


tol] 
bo] 


(xtax + ulguyat (3.2) 


Where the Se iia Mic moet bivc Semlacrinite matrix and 
RPE eeroOsttive détinite Matrix. An appropriate choice of 
mmeee Matrices must be made to obtain acceptable level of 
meee), «CC A(t) ©6=«60oand) «6(UU(t) 

The guidance law determined in Reference 8 is of the 


mocm : 


meen the performance index chosen (the weighting on the 


fergie vector was not considered in Reference 8. 





ee * ee ee 


on 
Il 
IQ] Re 
— 
wr 
f 
4 
re 
{ 
i: 


feteme tie tirst term is one-half the required miss distance, 
meeesecoOnd term 1s included in order to eliminate infinite 
Pees Olutieons. 6 as a2 by 2 positive definite, 


mame cric weighting matrix and of the form 


Bl 0 
B = (3.4) 
0 B2 
with ae aoe. wemoxImums acceptable value of 
- 
|U,(T) | ~ . A detailed derivative of the optimal guidance 


law is given in Reference 10. 
The solutions of the above equation from Reference 8 
Meye the optimal acceleration command at time T between 


merand Ti as 


Where ee Peewee immromrime tO 90 Util intercept and 
© 


a eweenow e-miss along the missile normal force axis and 
a 


moemiegative Sign 1S due to definition of the normal 


Gr 
On) 


MeeGeiwenaliOn a5 pOSitive in the upward direction. The 
epremal ,roli@rate command at time T between TO and fi 


meomeneterence § 1S given by 


2 
. 
= oan 66 (MJ 2 
Po = a (3.6) 
a” oT ae Oe 
c go 2 
finere Dies Peeve miss alone tie wine axis. Defining 
mee roll error: 
M 
be = tan *| Pe] (3.7) 
e 
mumalternate form for Po can Memoptaimed by daviding 
mametton (5.6) by Equation (5.5), then 
ae? 3 z 
ce Lae ae ae By) 
«oe aT (A$) eas) 
an LE ioe 1S. 
ce go S 


Berecuse the solution of oe was based on an assumption of 
Memiemissile roll excusions, the tanAd term in 


meat ton (3.8) can be replaced by Ad to yield 
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immoral cOmerolléeér 15 implemented in the present work 

Bewemomeawations (5.5), (5.7), and (3.9). In computing tne 

Sontrol input, ee >» Y and Z-component of P2ZC-miss and 
© 


Siemerememes Of the weighting matrix are needed. 


Meee APPROAIMATING TIME-TO-GO 

Ponte cescary £Or the Time-lo-Go and the Projected- 
Simos -)iStance in missile body axes to be specified at time 
Meo cCOMpute the Optimum control input. The Time-To-Go 
foe be defined at first, then the state variables will be 
Mmieroduced a5 means fOr formulatine the necessary equations. 
Meemmearoet acceleration components of the state variables 
meme be modeled Since it iS assumed that any information 
Meeut the target acceleration cannot be obtained from the 
active seeker of the missile. Then using this target 
Meeeweration model the algorithm for computing tne other 


Meate Variables will be developed. 


iiemlime-lo-Go until intercept, T; Ptoenedutreda TOT 
Mememtine the comtrol input. Let the time to intercept be 
Memined as the time to minimum range. Under the assumption 


that the relative velocity vectors may be considered 
mealteaie at 1tS present value until intercept. T can be 


mod at time, T=, as: 





pa a Ca) 





femens Rk 1S tne present position of the target relative to 


me missile and the factor in brackets in the component of 





Meeeone Rk . An alternate form of T 1s: 
es Poke ee = Roy 
T. = Senrek ory ia es 
1 2 7. Zz Cm 
1 \ a ee oy) 
io ry eZ 


The negative sign comes from the sign of be as shown in 


erire 2.2 and Equation (2.6). 


fee Oo ec PE W-ZERO-CONTROL MISS (PZ2C MISS) DISTANCE 

item) ro;eceted-Zero-Control miss distance (PiC-miss) is 
Bomemred as the miss distance which would occur if there 
MemeceenoO further control input after time T 

[iieeiis wOrk the State vectors are defined as the 


following state variables 


X = (R R R V V v a a a 
_ ( a ery rs” tx” : 2) 


Where R. Peeovemmcatiye missile distance to target 
ne. Pieomcc rot Ve Missile yelocity to target 
a, Pometcmeadroct deceleration 
Memodel £or target acceleration has been assumed. This 
Memciven as follows: 


a, (T) = a > 0 (Gor) 





Mirere a 15 the intial target acceleration and % is the 
Mmeemeprecal Of the acceleration time constant. The initial 
barget acceleration, 49's, will be asSumed for each 
Meemario, and « = 1/20 for evasive maneuver as indicated 
mn Reference 10 will be chosen. These values will be used 
memede Simulation of the optimal controller. 

[eeehethe assumption that a state estimator is available 
Mmempnrovide estimates of the state vector at current time (t), 
Mmmm aities Of the state variables for no control inputs at 


mere time, can be obtained easily using the integration, 


me. , 
te 
(ye Ga “J a, (t)dt 
iE 
R(t) = R(t0) +f V(r) at 
t0 


Meeeresuits of these integrations are 


20a 


we, ; : t0) eo 
Pees (tO) + 1/2 (1 - e vac d) e521 5) 


esimemccO) + (¢ = t0)*V (to) 
(3.14) 


: [1/o*(x(t Sele onl - “a(t - a, 


Meyere tQ and t denote current time and future time. ie 


Should be noticed that above integration is based on the 


Ol 





Peep eronm OL a Single physical dimension. If this assunp- 
mois 9s Foc applied to three-dimensional problem, it must 
Pemecteluidea that there are no coupling among the X, Y, 2 

finension. Let ec ee Epemecne COMpenents of Ehe 
Mea@omiss if inértial axes at current time (t0) can be 


eeepressed as 


M R ee 
; ‘ Tae 
M = R oF 
- go ey 
\ 5 
M tat Ro TO ee TO 
” 0 Q dg ] 
if 0 f, a (orc Sy) 
: Ey tte TO 
0 
: 0 
ae 2 
C2) ene 
2 Z z Oe -aT 
Pimere ¢ = f = £. = l/a’ [aT aol <e CC jee onedmtne Last 
a y Z gO 


memmeadccounts for the free fall effect of gravity in the 
Mertical axis. 

The components of the PZC-miss in missile body axes for 
memerol input calculation are obtained via the Euler trans- 


jmerrmation matrix: 


oS 





\f 


bx ex ps 

= ia | ra 7 \{ = |, M det 
se Lo} fe Tie] {A] ae ioe 6) 
lh [ M. no 


feeeecee tt 15 assumed that the missile has exact knowledge of 
its orientation angles vy ; a. Meee Communic accurate 
Mmeertial measurement units. 

ites valves Gt Optimal control vectors (a. and a ; 
men, are calculated from time to go cr.) eee hiss fain 

S 

omy axeS and weighting matrix components (Bl and B2) which 
[emebe determined according to the desired miss distance and 
M~e~meccootime time. in a later section, the weighting 


Beaeeex will be developed and the scenarios will be set up 


mom computer simulation. 


meee MULATION RESULTS 

Miewcuidance law in this work was tested on a simple 
memeeececorce-ot-freedom Simulation using only plausible 
mmm@eers without reference to any actual missile or design. 

Before the computer simulations, the weighting matrix 
mmm DeTrOrmMance index selected in Reference 8 will be 
Semeemraced Wwsing the initial state, X(0) , the missile 
Mererating time imposed, T (5 sec) and the constraint input 
Memmors. iff 1S assumed that the missile has a velocity 


@dvantage of two over that of the target, a zero- and 0.5 
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Beeemoe ime lace auto pilots for pitch with limits of -3s's, 
meee coe and fOr roll rate with limits of +6.28 rad/sec. 
Stemearrovoole miss distance which is usually determined 
Dy Sizing of the warhead was chosen as 5 meters. Which 
Meees the first term of the performance index Equation 
Meme caqmal to 12.5 meters. The acceleration term and the 


fewer rate term in the performance index are set equal to 


fms respectively (one-half of square mean miss distance), 


then 


2 é 
— s aa” ac = | 2355) m - (52.05) 


feenoninal time of flicht for the terminal phase is set 


Bema tO 5 Seconds. The components of the weighting 
Memenix B were based on the following assumptions. it 1s 
2 
assumed for (4g's) as a mean-square value of (a. - a 4) and 


(2 rad/sec)* meme an-sqtare roll rate to be acceptable. 


Mmieme the components of weighting matrix are: 


ve) 
ll 


Woe 254 cec” 


td 
it 
}—! 
ie) 
cr 
i) 


iis Ce 


tN 
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Meweolmmamas a and P are computed every 0.05 sec 
memarest1Qnen when the Time-To-Go falls below this interval. 

Mimotaet eso simulate the Optimal controller appropriate 
meemarios were Gevelopead below with restriction that the 
Meene=1O-GO WaS 5 Sec and the speed of missile is two times 
Mepeen tian that Of target. Three scenarios; tail-chase, 
Mee on, Side-approach engagements, were chosen as possible 
Maoes in this work. 

mets Meme tall-ecilase enIacement, the target flies 100m 
Sm@eve the missile and 2500m uprance with a horizontal east- 
Reemmeeceelerarion of 4¢'S. Figure 53.3 shows the case Z. 
Mmearse 5, Side-approach engagement, the target flies 400m 
above the missile and 1900m uprange with a horizontal north- 
mera acceleration of 4¢"s at t = 0 as shown in Figure 3.4. 
Memesimplicity, when the scenarios were set up, 1000m/s was 
taken aS a missile speed, 500m/s as a target speed in case 
1 and case 2, 600m/s as a missile speed, 500m/s as a target 
meted in case 3, and 4¢'s as a target acceleration for all 
three cases. When simulation was executing, this first 
MeteraetoOn assumed perfect state feedback without the use of 
emeootimal estimator. 

Mirensimulation results are shown in Figure 3.5 through 
Mere 5.10. At first, the optimal controller with no time 
Memedomitested fOr three scenarios, as shown in Figures 5.5 
@meaeo.0. For scenario 1, the normal acceleration command 


mimemictial phase is very large, then decreases to zero at 
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remo. poesia tO>] TAIL=-CHASE ENGAGEMENT. 
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Sone eee NNO es HERB ON ENGAGEMENT. 
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Mmeourkc 5.4 SCENARIO-5 SIDE-APPROACH ENGAGEMENT. 
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ommelle) find! time. For scenario 2, the acceleration 1s zero 
Mumem@oniNnitial phase, then increases and converges to zero 
mmeeiy. FOr scenario 3, the acceleration starts from zero, 
Mieeiicreases rapidly up to the maximum acceleration tnen 
Mmeeees £Or Seconds before decreasing. The acceleration 
MeeemeiOu converges to zero at the final time. Comparing 
Mebefe  GeSULTS With the variation of the roll rate command 
Mempeboure 5.0, Lhe Opposite trends are observed in the 
jmemcs Of scenarios 1 and 2, 1.e., if an acceleration 
@emmand in the initial phase is large, a roll rate command 
meomeal!. tlhe acceleration command in Equation (5.5) is a 
Mme uGerlon ©f PZC-miss in of z-direction of the missile 
Mmeerdinate System. For scenario 1, PAC-miss in the 
mmemrection 15 relatively large, this causes a large accel- 


Mmeteron cCOmmandgd and a small roll rate command in initial 


iT eeferrsceiario 2ytne relative magnitude of the PZC- 
Meme che Z-direction is small. This cause a large roll 
meeeecommand in the intial phase. For scenario 3, the 


Mmeemiss im y-direction is of the same magnitude as the 
Memative distance as the Time-To-Go decreases. This means 
that a larger acceleration command than the maximum value 
femmewablie in the missile body coordinate system 1s required 
Memeeome Of instance time during the missile flying. Thus 
m@emeccelcration does not converge to zero at the final time. 
in the derivation of the equations for control commands, a 


small angle approximation was assumed. However, in the 
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eee o 5 this assumption is questionable. The roll rate 
command computed as shown in Figure 5.6 shows a sine wave 
Garacteristics which is possibly due to overcorrection. 

Ghe simulation results for a 0.5 second time lage are 
Sugg 11. Figure 3.8 and Fisure 3.9. The time lag causes 
Mmem larger acceleration control commands. This includes 
Meercorrections, so that the Eencrol commands “do not 
Meee ree in the final phase. The other characteristics 
memeeow the explanations given in the analysis of the cases 
with no time lag. 

Pe colietaree ScCelarios the miss distances were obtained 
Mmoecide Simulation of both cases, no- and 0.5 second-time 
meee ihe miss distances are below 0.5m for all three 
feemiarros in the case of no-lag. In the no lag case the 
Meer ol laws gave excellent results. In the lag case, the 
Mumemectoeances are 5.9m for tail-chase engagement, 0.8m for 
free on ehncagcement, and 32.9m for side-approach engagement. 
Mteemiss distance of the side-approach engagement with a 
Seeeeerecnt maximum acceleration limit fo 23¢'s results in a 
Memeemecistance of less than 5m. This higher maximum acceler- 
meom Jimit will be used in the scenario 3 during the simula- 
memes fOr testing the performances of whole control system 
Me cnapter 5. 

W@enorevyious results are summarized as follows: (1) the 
im@tenine Of a target is very dependent of the missile and 


Mmemmerety peometry, (2) the capabilities of the missile 
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wiemum accelecartion and maximum roll rate limits and time 
Bolstant) in maneuvering against the target is very important 
Memeeers, (5) hicher maximum accelerations are required in 
Mieside approach engagement, and (4) the assumption of 
mel Ad in the optimal control law is a limitation in the 
Memeral Case. Nevertheless, the control laws developed 

mom the Optimal control ™theory for both the no lag and the 
mmmmeerecond time lagucases yieided successful results for 

Memeo gect7icd SCenarios. Since the optimal controller with 
MepemcccOid- time lace auto pilots resulted in amiss distance 
mimeoproximately Sm with complete state variables feedback, 
MemeeeecxDectéed that the inciusion of an optimal estimator for 


Meteniarricular scenarigs will results in still larger miss 


mrs tances. 
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fee DESCRIPTION 

Teenoue tui) State variable feedback it 1S necessary to 
eeemace the state variables. This chapter describes the 
Mmeman tlilters used in the state estimation. These filters 
Memeemcecsicned for respectively: (a}) two measured states, 
Beye tOlir measured states, and (c) six measured states. 

Ali measurements have noise superimposed on them. It is 
@eoumed that the nature of the noise source is Known 


Semoletely. 


imeaijameriiter 1s the optimal estimator of the 


[eecs. olnce the measurements are non-linear, it was 
Memeesosary, tO use the extended Kalman filter theory. A 
meee Orcer extended filter was assumed in this work. In 


See lOnine the filter algorithm, the non-linear measurement 
eeiteton 1S linearized about the most recent state estimate 
fmeeecien the xnalman filter alsorithm is applied. This is a 
Seeeemded Kalman filter algorithm. The results of the 
meeoded khalman filter will yield suboptimal state estima- 
iron) . 

Mempencetiece, elve missile will have sensors of various 
i[epes that provide inputs which are related in some fashion 
jmomeie States to be estimated. The main emphasis in the 


DBresent work will be on active seeker which is assumed to 


54 





provide the optimal estimator with measurements of line-of- 
meet ancgies for the case of the missile ay Wie we ite as re 
Meme ocnsors. For the case of four measurements one has two 
of line-of-sight angles, range and time rate change of 
meme. FOr the case of sim-measurements one has two of 
Mmeeeot Sicht angles, two of time rate change of the line- 
Meee cgit angles, range and itS time rate change. When the 
feeetirement equations and the Jacobean matrix are developed 
meee x Measurement case will be vderived since the other 
Meo have the same Kind but less number of measurement 

ems OrTs. 

Peele first consider the state Equations and then the 
Mereeurement €quations. This will be followed by a discus- 
memeor the extended khaiman filter. Measurement error 
Memeeets Of the initialization of the kalman filter will be 
femezed. Finally, the estimator for the measurement cases 
fameeoe Simulated in order to test the effects of the imple- 
Memeation of the different measurment vectors on the 


extended Kalman filter. 


fee oO LATE EQUATIONS 
mieestace vector tS limited for convenience to the 


following state variables: 


: a 2 
Pa 


ie 
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Phere XV Rv, A, = R', A. = R, are the components o: 
“ — ae - a 
Meme vye position, \ = V__, +A. = V.., xX, = V tne com- 
+ Ee 5 ry 6 re 
Mee ges Of relative velocity, and xX = a,., \, = a._, 
/ tx § ty 
Ng me oe, HicmeGupolicmts Or Carget acceleration. 


iMumm@rcetT LO Cevelop tie State equations, dynamic equa- 
mmems OL target motion are needed. The target model 
Meeeected for tracking applications must be sufficiently 
ie LO permit ready implementation in weapons system for 
Meme COmputation time 1S at a premium yet sufficiently 
Sememlsticated to provide satisfactory tracking accuracy. 
memeet tiese requirements, the model described by Singer 
Meet. 10] wiil be used in this work. The model will be 
Mm@eeented fOr a single spatial dimension in order to 
mele accurate tracking performance estimates to be made 
mma Variety of sensor measurement. If the targets under 
Semerderation normally move at constant velocity, the 
Meeelcrations due to turns or evasive maneuvers may be 
Mrewed aS perturbations upon the constant velocity trajec- 
Meee ihe target acceleration a(t) therefore wili be 
memmed the target maneuver variable, in a single physical 
Me@etston. the target maneuver capability can be satis- 
Mmemeorliy Specified by two quantities: the variance or 
Mememicuce or the target maneuver and the time constant, 
Memmeotmration, of the target maneuver. Hence the target 


meeerctration, namely, the target maneuver is correlated 
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Mempetinc- A typical representative model of the tarc®t 
Meee retauign as a2 correlation function, can be expressed 


we, tOlliows: 


2,7a/t| 


eet) 5 Ela, (tja,(t + t)] = ™ a > 0 (i) 


nN) 


Migere the On Vem itawecmot the target acceleration 


mame 6ciS «6the reciprocal of the maneuver time constant. 
? 


— 


ime variance oo Gusetenmodelewill be approximated using 


meer ollowing formula which represents the acceleration 


Meee abdility density suggested by Singer |[Ref. 10]: 


Z 
2 omax 
GO = = ) ere aoe Cacd 
m 5 max 9 : ) 
where a Pesca mMawimim acceleration rate, P wees. the 
max max 
eeopadility of maneuvering at ais and Po iS elec 


Meeeoebiiity of no target maneuvering. And as mention in 
M@eeerous Chapter, a = 1/20 for an evasive maneuver will 
be used in this work. 

P@e abOve process can be modeled as the output of a low 
Mees fliter driven by white noise and it can be described 


ieee ne equation 


ae = - a. + W (4.4) 


oy 





Metere= co 6(t) «6 «vuthe: Correlation function of the white noise 
WV 
miout Satisfies 


* 


oo Go. Zag, 8(1) (25) 


W 


ireeaecceteration model above will be applied to each 
feo, 2 dimension with no cross coupling, the system state 


Semacions can be written then as: 


eee foe 





R. 00010 00 0 ojfR | 0 0 | 
XY | x | 
R. 0 0 0 0 1 0 0 0 OFIR 0 0 
yy 
R_ io og Lo 0 0 R. 0 0 | 
V a OCU re ODay -a | 0 
es Hee ay 

N woo DUC COCONC OC CO. Sica Baz 288 
ry ry | 

\ imo. 0 8 0 0 O TilV.. a. 0 
eZ, eZ Ve 

| 

ao . ee OO 0 2 0 Oita. 0 iW 

lw 

-, io 80 0 CO 0 0 : Olja,, | 0 WH 
a. ieee oO 07D 0 Be ae | 0 ey 


(4.7) 


fmemecquation (4.7) is a linear state equation with the 
mummmuron that the forcing function a vector, 1.e., missile 
Meeeretation, can be precisely measured and also resolved 


meee NX, Y, Z components in the inertial frame. 


| 0 


Mam ~eoensOrs ave a constant data rate, sampling the 
mera ev,ery | seconds. If the above Seciwecuat 1 On yi 
Mmeeresented in discrete equations. The appropriate system 


Seacions of motion may be expressed by 
X(K + 1) = @(T)XCK) + BC(K) + UCK) (4.8) 


mmete (1) is the target state transition matrix, B{(h) 
Meee terministic forcing input vector which is composed of 
mist and second integration and U(k) the inhomogeneous 
Birvincg input due to white noise. Since it is assumed that 
mmemiissiie acceleration is a known deterministic forcing 
mimetion, the system equation of motion can be obtained by 
@irect integration in each single dimension with E|]W(t)]|=0 


Mgeswaesired integration in discrete form yields 


R(t0+T) 2 Vy enac- 1 + aT + e @') | IR¢tO) 
V(t0+T)| = jo 1 a es V(t0) 
an(t0+T) 0 0 a feet] 
SS — Peco ciedt 
Cue | 
+ i an ee che 
0 
puke etl) =o emo eth) ee a (xk) (4-9) 


Where K+15 (K+1)T=tetT 
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, K = KT = t0 








Sete Pt tne drlivine input U{ct) vector is considered’, 
mms input if not a sampled Veroeonwee mr etme sCONnLINUCis time 
Meee mO1Se input WN{(t}) . Since W(t) is white noise, 
fie aiso should be a discrete time white noise sequence, 
Meets, 6&1 U(K)U(K+i})} = 0 for i 70 . Then as shown 
Memeeercrence 10, U and W are related bv 


a) 


moe Oo) iki Tt - 1] W(t)dt (4.10) 


kK 
fete U(k) is a white noise sequence with covariance 


matrix 
Couto ae (4.11) 


Pomlowine the abowe derivation, it can be easily 
Meeeemaed tO the case of three independent dimensions in Xx , 
Mme rOor O@ -@ =a. =a , and o_ =Fc, Fo, Fa 

x v 2 ae Vv z 


Pease tie model of target acceleration is for one dimension, 


mime transition matrix 1s 


¢ =| 0 I Pee (4.12) 


60 





Meeyvows. and 0 are 5 by 3 identify and null matrices, and 


> ae 
tr = Was (ct ec ae 
mee yeti - ¢ **) 
ee = a al 
2 


Meeeerorcine acceleration vector is 


aaieeee t) 
B e J - 
il KT ae or) dtdt 
Cheec weal | 
B see eet ace dt 
Z eT my 
loge ae 1 
Be ec de eG Ghe 
= KT Mz 
: ome wot cae 
B = ara a a ele 
- KT mx Gos 
(Cee ta 
Be er omy (t) dt 
ie oa) | 
B - cent © 
6 on Mz 
BK 0 
Bg 0 
By 0 
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The covariance matrix of white noise U(k) is 


Qik Ql st | 


g - | Qo! Q551 oot (4.14) 


Where 
1 = Cor/a") tL - e888 + 2aT + (2/3) (aT)? = 207)" 
are 2 | 
Qo = (o7/a>) fe *F + 1 - 2e°%T + 2aTe 8? - 2aT + (aT)*] 
Sf (37 /a*) [1 | ere sa 
Q5> = jie 7 ore! 2a] 
Q55 = (s2/a) [e287 + - 20787} 
eR = aie 7 eee 


eee MEASUREMENT EQUATION 
Pmaomed>e OL tic CStimator with six-measurement sen- 
femeemeene present study assumes that the active seekers of 


meee Giye the measurements of the sightline angles 





Sp» Vp > Sten Vesta wcChange Of the Sight Jine anylés 


8p , Pp PovUCmuciativesmainoac Rk | and the time rate chance 
2 ‘ 
Gf the relative range R contaminated by the Gaussian white 


Meaoe. Then the equations of measurements in discrete time 


form are 
Z(K) = h(x,K) + VK) ers 


Merere (hk) 1S a sequence of measurement vector contaminated 


Dew hite noise V(kK) and 


| 
ete vate ame e = | 


r 
Seemateren wanele rate ee 
azimuth angle re 

h(x) = = wa 
azimuth angle rate i! 

= } 

| 

range R | 
range rate R | 


Mmeieosile seekers actually measure the data about the 
Meemey axes. Thus, throughout this development the assump- 
Meee, iS made that the missile possesses an inertial measure- 
Mamemunitc that accurately specifies the missile's orienta- 
Meemein Space so that a transformation from seeker to 
Mmm@ertial coordinates can be made. 

Deterrimg £O Figures 2.2 and 2.3, the six-measurement 


Meerors can be written 








ae ie Ea aS Ar ene aaa 


wy 
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but 


eZ 
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a an ~ 
c ey eae Oh yy a) 
(ee 
x ve 
a rRV +RRYV 2 = pe 
! - 
Woe ak 2 Ys y Z ry ( x y )\ 
R 2 2 2 2 Doe 
(R + R + R BGs “se RY ) 
R 
-] \ 
tan 7 
~ 
ae ee +e RV. 
rw z y 1x a 
a 2 ‘ 
R COS¢,, ie = RB 2) 
x x: 
Z ? 2 
CR. + Ss = Ls ) 
ee Ne 
B-W/R - A a 
eee oe ieee 
ne y z 
Ro = Ny ; ae X, : R. = Xs A one = Ny 
= Xe 





imewmedsurement vector can be expressed in terms of state 


mectcor 
x 
-| me 
phy | re 
| | 
| oe eX. + [x ex 
Ih pee 2 ss aa 2 2236 | 
2 eri a. see. ay | 
Pia + xy — | | 
i 34, 
h(x) = ia. = tan — 
—s— S x 
1 
| one 
4 er meaner | / >. 2 Re 
(Ay = = a As ) (A, = ‘5 ) 
; - | 
| 
Dea 2 meer 
he CX, + XN, + \ , 
SC a ae | 
h _ 
eo een) 2 
| % | (X,° + X57 + x,4)7/ | 
(4.16) 


ites measuremene moise Covariance matrix for the active 


seeker can be written 


2 
Tp 0 0 0 0 0 | 
conn ce 8 «(0 

" eee 0 0 0 

E[V(K)V(K) ] = R = p 5 (4.17) 
0 0 0 oy" 0 0 
0 0 0 0 cae 0 | 

sf 
0 0 0 0 0 on | 





mmetewtic diagonal elements are the variances of the indi- 
mecual measured quantities. 

As shown in ieee oom che measirement vectors are 
Memon -tinear f£Mnctions of state vectors. It is impractical 
memrmplement in non-linear form because the computation of 
Meee KR and GCrror CcOVariance matrix P in update process 
[Mieextended khalman filter algorithm is not possible. To 
emery this Computation, and to implement the extended 
Mmeemieritter, the Jacobean Or matrix Of partial derivatives 


feewill be determined, 1.e 


- 3 


Hee > ay (4.18) 


mmienvacobean will be a 6 by 9 matrix because h is a vector 
Memesix ciements and X 1S a vector with nine elements. 
Meeetorming the operation, the components of the matrix H 


ear C 


H = ~X XK Ye ton ee) 

Meee gs 5 

5 2 

EH = {X Xx ite 

ud \ 71 gi7 via 

E = 0 4H = Omen = 0 4H = 0 
U6 u7 ué 4g 
H = (X )/(D 

64 ( A ( c? 

5 =X » 

Wee ( 5I7 | ) 

H = (X D 

£3 ( 3/7 c? 


66 





tl 


0 4 mere llCUCElU(DCUKC OD as 
3 Sb i 8 rod 
? Q 
iets x) 1700 
Bone e ee tess 7D) 
lee : BD) 
xX = ox U 
aa 7 5 ma 
2 , 3 
. ae 
[OO %) - Oe) q7 8, ) 
ao) 
1 - 
x > 
es ! 
ony (DD 
3 0? 
ioc eeeno | 8S = 9 
68 69 
Z Z Zo NZ 
(X Pexore ek.) 
Z a 
2 ge 
(A PhS) 
ie 
a a. 4 2 
ex [cx Pano Oy =< Xe (x +X )] (30 +D 
He eG 6. 4 1 2 SG 1 0 


ek ( x ; aa es D : X : x iy +k eg et os 
i +3 L [ = - Cc 
2 sl = 6 71 > Jd j C 


2 2 2 on 2 
ne [XAceS ok + xX Xx + xX D +? 
(x, eae ms a] ®, soe 


Z 
oy oD) 
( alles ; ) 


? 


Xx X 1D Oo 
rae 0 i) 


D/L “) 
1 C 


Oe? Wes = 0. 4 = 0 H = 9 H = 0 4H = 
hd 16 17 eS 9 


67 


2 


) 





tr) 


tn 


if 


i 


ad 
bra 


( 57 | C 7 


peed. eA, 
172 3 4 


uy 
) 


Nc (D *5 : 
x pay 
5 J Oo 7 


68 





Meee the EXTENDED KALMAN FILTER 

eae eine eOrsmelacemiiimize the estimation error, 
Beamebe divided into two processes, the one is the updating 
Mieco TO cStimate the State vector at the current time, 
mased upon all past measurements, the other is the extra- 
Mepereinc process to estimate the state at some future time. 
Meeinal estimation procedure for linear state equations and 
Seeeinear Measurement in the form of discrete time will be 
Mepreribed yererring to Reference 6. Fisure 4.1, a timing 
Meeeeram shows the flow of the various quantities involved 


Mmemiene discrete optimal filter equations. 


1 


Meewaitscrete System Equation whose State at time kt is 
Memotcd by Simply X{k) , where U(k) is a zero mean, 
Mepeee Sequence of covariance Q(K) , is 
PC oe BR 1) +t UCR = 1) ia) 
Mrere © is a transition matrix 
EU Ce =o ee 2) 
: a a 
Ee UC iC 5 ) J = eG (4 3.25) 


mime measurements are the non-linear function of the system 
Memes Variables, Corrupted by uncorrelated white noise of 


covariance R(K) . The measurement equation is written as 


69 


ZK) = h (X,K) + ¥CK) #24) 


where ELUC(K)] = 0 (1,25) 


ECVG)YG)'I = R8,, (4.26) 


mime initial conditions are 


> ee A 
Meco) ] = X(0) E[(X(0) - X,)(X(0) - %)] = P, ee 7) 
aes | 
EGY Vy ) = 0 uncorrelated (eee 


iewecilean filter alcorithm is tO minimize the estima- 


as . 
Mee Ot emis c., ciror Covariance. If X denotes an esti- 


dat 


Fictional 
Memeo. State vector, i , and X is the mean value of the 


cr 


mmeee Vector, the error covariance matrix is defined as 


Bee OT: . 
Peer eee C= NX) (symmetric) Ga. 29) 


_—— 


mmemdevelOpment of the extended halman filter algorithm with 
epee detinitions 1s essentially more application of the 
expectation operator Pn eeva Mieonmecider the tpdating 
Mmieess then the extrapolating process. The update equa- 
mmeons are used to incorprate the latest measurement in the 
ememilace and in the covariance. After we obtain the updated 
Mmeariance matrix we will then consider the optimum choice 
Seekalman gain and derive the equation for the esti- 


Meeweommor State. Let X{k-1)(-) denote an estimate of X 
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mone) R{h- 1) 


X(K-1)(-)| XQ-DC) 


k= 2) ,O(kK-2 Oe serra!) 


> 


Mee - 1) (-)} PG. (=) 
erenameka tlie 
process 
ers) 
measurement 
pueeecss 
MeeGURE 4,] 


al 


V(KYE -) le) =} 
eat } 


2{h),Qts) 


ae 


el) ees) 


¥ , 3 fo ‘OF 
UWDac cing 


process 


t (Kk) 
measurement 


process 


DISCRETE EXTENDED KALMAN FILTER TIMING DIAGRAM. 





eee wn - i 5 Given Measurements up to and including tine 
Mem (h-1){-)  eenote the error covariance, K(k-1) 
Semote xnalman gain and H(k-1) denote the Jaccobean matrix 
Meetime K-1 . We will seek an optimal observer in the 
Meeesoence Of the state excitation noise and in théepresence 


mecne opservation noise. 


H 


P(K-1) (+) E4C1 meee Me i pe CX ¢K=1)!) ¢-) (1-K 


Pee eh) KCke2) | 7 


+ 


emer unoG oI ih-1) (-) | (1-K(K-1)H(K=1) 2) 


+ 


v(k-1) 'k¢x-1) 7a} 


Peteeamerme Lauacmon (4.50) using the definitions 
Meme ions (4.26) and (4.29) and the assumption that the 
Measurement errors are uncorrelated, then one has 
fmix-1) (+) = memento 1) 12( k= 1) (-) (1-K(K-1)H(K-1) |] 


(#31) 
P@@e) RU 1) K (K-21) | 


Awe ouartOmeupadtes the error covariance matrix. 
miemewe Still need an Optimal gain kK. We will choose to 
feemize the diagonal elements of the covariance matrix, 


nee. , 


J(K-1) = Paar) (2) CK) (#) 1 ate e(koty(+)) (4.52) 





Mamememnmatprcadi- mode | 


—e oe eee ee ee ee EE ee eee eee eee eee oe Geet ees eee ee ee Ee ee ee Cee 


mm rn re mm ee ee ee ee 
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MeeuRE 4.2 SYSTEs MODEL AND DISCRETE KALMAN FILTER. 


re 





Piero Veimie sd Dit Of Matrix calculus, one has 


Cy] 
O70) 
Nee 


5 berace A B ay 


2AB ide 





se 


[-K(K-1)H(K-1)] = -H(K-1)' es 


ae eee [ 
ak (k-1) 


Ppplving these formulas to Equation {+.5i} and solving 


mor kK(k-l) , 


Mey | PCK-1) (-)H(R- 1) te) Pe 1) ( -) H(K-1) '+R(K-1) 77} 


eS) 
mes 15 the Kalman Gain Matrix. 
Pearranging Equation (4.51) 
P(K-1) (+) = [I-K(K-1)HCK-1) JPCK-1) (-1) (4936) 


fiiem the update State eStimate at time K-1l is Equal to 
Meee xtrapolated state estimate at time K-2 plus a term 


which weights the measurement residual via gain K(k-1) 


A a D> 
Meee) =~ NtR-1)(-) + K(R-1) LE(K-1)-h(k-1) ,X(K-1)(-))] 
etme 7 ) 
Bomiced tO CONSIdcCYr also the propagation of the estimate 
es between measurements and the propagation of the 
COVariance matrix P({k-1)(+) between measurements. Propa- 


Meeron Of the estimate is straight forward and involves only 
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Ee eanG ) 
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PStimatre 


eles A 
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GOmMpute Mat) mea cSomanlenesss ie 
ae 7 q Cnecks 






n(K) 2 compute X,(-) : 
4 Peet) 1 


TE I rE BS eee eee oe Fe ea es 


Meek 4.5 DISCRETE KALMAN FILTER INFORMATION FLOW DIAGRAN. 


is 





fhe application of the state transition matrix and the 


mereing terms. 


ws a : 
RERIC-) = O(K-DYRK-1) (4) + UCK-1) a 


Cs} 
CO 
Nee 


Peer wcenmnuenem (29), solving for P(k)(-) 
; a : ; il : ae 
Poe ol PY Re) te oCkh-1) «+ QCK-1) Veo 


Pu@iiGema Thidwseradtes the above developed equations in 
miock diacram form which contains a system model, measure- 
Meme process and the kalman filter to be implemented, and 
mere 4.5 Shows a simplified computer flow diagram of the 
meescyete halman filter. 

emer eos eelroOuch (1.59) constitute the recursive 
formulas for implementing the extended Kalman filter 
Meeorithm. the process is initializing by providing values 
fn 
Meee) and «p(0)(-) . In order for the extended Kalman 
meltcer to be simulated, the next section will describe the 


Memmeetralization of the Kalman filter. 


Meee INITIALIZATION OF THE KALMAN FILTER 

i womemermcerlyvewOrernis wOrk 15 to test the perform- 
mmeenOt Bli missile control system having a kalman filter as 
fmeoptimal eStimator as a function of different measurement 
Meemors. All values for initialization of former worker in 


Meese work (Ret. 14) will be used for computer simulation. 
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Te 


eee iene tie error Covariance Matrix, P 





a 


bet a, » Cy and oo @e27 fe Che Standard deviatyon 
Mumeeeative position estimate, of relative velocity esti- 
Memeo and Of target acceleration estimate. The initial 
Memees Of Error Covariance matrix p , will be chosen as 


mmerows. the first six elements of the initial state 


Meetor, are determined by 
Nn 
X(0) = X(0) +# V(0) 


Metre the measurement vector x Peeves t nucemeas Urement 
meeeor \ plus ¥ caused by measurement uncertainties. 
Smee xpected mean square of initial covariance matrix 
Meoonenits. with the assumption that the initial 
Memeariance matrix is diagonal with the diagonal elements 
Memeeecting the uncertainties associated with the initial 
MepmeemvectOr estimates, the components of the covariance 


jeer ax are 


2 ees 
oy te - 2%) = EY, ] =, 
r a y z r Z 2 


Mmemerast three elements of the state vector, acceleration 


MemmpenentS, are set to zeros 


Kg = X, = O 


tod 
l 





Pen letalizine Une Mreaslirement Noise Wariance 


—_——_— > er a Sed —— ee 





Measurement errors tor the active seeker are 
meoumed to be due to thermal noise, gimbal angle pickorr 
omeor, environmental noise, and glint. The values of one 


meme eyrors ihn each measurement will be used in present 


fi 


Mere because in the G@ase of angle measurement, the envi- 
memmental noise iS proportional to the square of range to 
imeeeet and glint errors the wander of the apparent target 
memmeroid as seen by the seeker as varing inversely with 
Meme tO target. Ihe errors Caused Dy noise are generally 
mmeenl, complex function of the target geometry, target 
Meme CTOSS SECtion, radar receiver, signal-to-noise ratio 
emeeto. FOr Simplicity, in present study, the constant one 
Pm errors (average values) in each measured vector are 
Mewmed tO initialize the Kalman filter. Referring to 
Meeeerence 15, the one sigma values of the average errors 
mememeaSurement Of angle, angle rate, range and range rate 


Meemene present active seeker are tabulated below. 


mel angle measurement errors Peek Smee On deo 
(2) angle rate measurement errors : 0.5 -2.0 deg/sec 
Wo) range measurement errors > meters 
(4) range rate measurement errors : 6 m/sec 


Pore soliibotton runs 1S Undertaken, it will be 
feeeer iy discussed the manner in which the Kalman filter on 


Memeo cine Missile is initialized at the start of the 
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engagement based on sensor data from the launch aircrart. 
BepmelommcmealtCta:t'’Seseisors are assumed capabie of dscter- 
Semato, Within sOme prescribed accuracy, the relative 
Meottion and relative velocity of the target with respect 
memene Missile at the instant of launch. This infcrmation 
ewe POmiiahiainee tie first Six €léments of the filter's 
Meeate vector. However, the launch aircraft is assumed to 

be able to provide no information regarding target 

Meee oul. lms the three elements of the initial state 
Memeo; are set to zero, finally, the complete list of 


mominal parameters for initialization is provided in Table 


ei 


See AL UATION OF THE OPTIMAL ESTIMATOR PERFORMANCE 

ites peritOrmance Of the state estimator with two- 
measurement vectors (case 1), with four-measurement vectors 
(case 2), and with six-measurement vectors (case 53), was 
merated tO evaluate the effect of the possible implementa- 
meom Of the more measiirement sensors on the bank-to-turn 
fmeorle £Or typical scenario (tail chase engagement case). 
Miemerror CcOVariance and the kalman gain components 
Mmmeecrcd, and the states estimated, were computed as shown 
Mmemeeroure 4.4 through Figure 4.30. The control laws imple- 
Mented with the estimator also were tested to evaluate the 
Meresotmance Of the control system of the BIT missile as 


meri in Figures 4.51 through 4.35. 
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foemeeeaee Variable. Figures 4.4 
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error covariances o 
Mmpoug) 4.7 show the three covariances for the relative 


Mmesitions, 1.e., Pay 5 ee Pe - Only the covariance 


G! 


Seethe two measurement case varies substantially. Thus 
@menwould expect that the gain components G14 ; G54 ; G4 
memid show Only variations in the two measurement case as 
Muemeeivowi in Fioures 4.13 through 4.15. The trend noted in 
Mmmewset Of graphs (4.4, 4.5, 4.5) and in the corresponding 
iemMan gains 1$ also observed in the other covariance 
Memeoments. thus the covariance components for the velocity 


i 
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mamemres +77, 4.8 and 4.9) should larger variation is 


Mepeeocting in the corresponding Kalman gains of 4-16, 4-17, 
Smee ine COvariances Of the acceleration components show 


Meeeiidar development but there appears to be less variation 


an aay and a acceleration covariances. In the case of 
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a. eeomveridtetoneor the halmanm gains are restricted to a 
fmen narrower range. The error covariance and the Kalman 
meen COMpOnents in all three cases eventually converge to 
Meeeeeecomall values. Figure 4.2% through Figure 4.30 show the 
[meets Of the state estimations of the three estimators. 
Mmmeeures 4.22, 4.23 and 4.24 show the variations of the 
estimated relative position Gee Rio» ei Figures 4.25, 


feo and 4.27, the variations of the estimated relative 


mewocity Oe ee ee eee sme So  e4ee29 and 4.50 the 





Variations of the estimated target acceleration Ce 


o> ieee) Peet Urewemrorewtne curves representiny the 
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Mctimated states are similar except for Figure 4.25 which 
Bows a wide variation in the estimated ve Someone uloe 
jeme two measurement case from the true state values. 
Memeraliy, the estimated states of the estimator with 
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mao- and four-measurement are close to the true states, 
mmmeopite Of the wide variation of @rror covariance compo- 
Memts Of twO-measurement case. The estimator with six- 
measurement vectors generates the estimated states which 
show almost same characteristics as the estimator with 
FOur-measurement vectors, but usually underestimated the 
Me@ees i11 COMparison with the generated states of the 
meemmacor With four-measurement vectors. Figure 4.32 and 
metros. 54 Show the curves of the control commands com- 
puted using the estimated states. The results of the six- 
measurement case have a wholly different form from the 
mieves representing the computed control commands in the 
Other cases. The normal acceleration commands computed 
Mmemomcic CSLimtead States with the six-measurement vectors 
initially have smaller values, but finally reacn maximum 
Meee 25 Shown in Figure 4.352. Also as shown in Figure 
Meee the roll rate commands run from the minimum limit to 


Sieemaxamum limit. The large variations of the both control 


memmands Of Case 5 during a short time interval result in 
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Peete toto tie Six-measurement sensors containwani'™ 3: 


Meee er eel Orsmeiicl are exactly same kinds and ngimeers 


ih 


mee Measuravic state vectors using the rour measureyvent 


MemsOrs in case 2. As the angle rates, 


YR 2 op 5 
mmermidaed to the estimator as the measurement vectors in 


Mees, the both angle rates are non-linear functions of 


i 


v ar es ar 
a\ aN and ax 
eee ae 6 


Mey SiX-State vectors, Ny > N55 As 
Me ociitS Of the G€stimator with four-measurement vectors 
Mimeewitil SiX-meaSurement vectors show the almost same 
Memematcions Of the error covariance and the Kalman gain com- 
Meer es Of the SiX measurement case as of the four measure- 
Memmencase aS shown in Figures 4.4 through 4.21. 

eomee hand fllter in tCmis work 1s a first order filter. 
Mepeeehe COmponents of the Jacobean matrix contain only the 
Mueeercer term, 1.¢6., all higher order terms are neglected. 
Memeo Order of the system is increased one would expect a 
Meeeemcompiex svstem to response differently. The above 
Mer-Orls attect to the underestimation of the states of six 
meeeeurement case aS shown in Figures 4.22 tnrough 4.50. 
Meemmetie above Considerations it can be concluded that as 
Mim@enwincreased number of measurement vectors is implemented 
Mmemeene Extended haiman filter algorithm, the result of the 


Meeeem May not be enhanced up to the expected degree due to 


Meemerereascd COMplexity of the system. The system with the 
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order filter as a optimal estim 
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epeeese NOW se ewcreeies the performance of tne control 
meecem implemented in a bank-to-turn missile for each scen- 
moro, tie maXimum normal acceleration, the maximum roll 
Meenrand the time constant are the key parameters which 
Meeerect the missile Capabiiities, the noises to measured 
Mereeors are the parameters of the environmental effect, and 
Miliary variables of some importance are then total en 
ment time and missile velocitv. The mean miss distance 
@etermined trom the 50 Monti Carlo runs was used as a 
Meee ormance Standard for the estimators of two- (case l), 
memi= (case 2) and $Six-measurement vectors (case 3) on each 
Meee oO. Ine Simulation results are shown in Figures 5.1 
mememen 5.19. The sensitivities to variations of the para- 
Meeers Will be discussed below. Since the simulation results 
mm@eweessentially the same characteristics for the tail-chase 
engagement (scenario 1) and for the head-on engagement 
(scenario 2). The performance on these two scenarios will 
Memeotscussed together for each parameter. The side- 
Bporoach engagement (scenario 3) will be discussed sepa- 
mately. 

Peeures 5.12 and 5.2 show the miss distance variation as 
Peeeeiyet1On Of missile maximum normal acceleration. For 


4 


feetario 1, case 1 shows a mean miss distance of 9.8m at 
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Bemere loss, Case 2.9m at above 1l5¢g's and case 35, 15 nm 
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Meeetne normal acceleration decrease 
meee is 1eS5s Sensitive than in the other to Cases. 
Peeescenario 2, the result of case 1 shows tne minimum mean 


Mm ons ance is lO0m and the result of case 2? 1s 17m at 


x Cees Fee pectivew a tite result of case 3 is 32m at 


Sax eee (emo wetoemom ther curves show that the mean 
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m1 Peeee nce 6 CasSCulwand Case 2 Wary in same manner, 


Nh 
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Meme rises Sharply as the maximum acceleration decreases. 
Seecne Other hand, the mean miss distances in case 3, vary 
smoothly above Nees Pie nay Conclude from the 
Mien sis as fOllows: for scenario 1, the mean miss distance 
mame decreased since more information on the state vectors 
mumeie Obtained from the increased measurement vectors in 
Seeenacor fOr low maximum acceleration. In order to get the 
Mem miss distance below 20m for all three cases, one needs 
os x ae ee cero SCemarl@me as the measurement state 
Meerors are increased in estimator, the mean miss distance 
Mepeeot decrease, possibly due to the effect of neglecting 
Meer eher order terms in the extended filter. As the mag- 
Meeeoe Of the values of the missile and target geometry 


MemeonentS becomes larger, the effect of the linearization 


becomes more pronounced. eee over 13g %s is required to 
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meres 5.5 4nd ov Ene Varlations or the mean miss 
MMPonGes 425 Missile maximum roll rate varies. Lie Seen aro 


Mumeror case | and case 2, a mean miss distance of below 10m 


Mme be Obtained when Le Mepemeiter than oS rad/sec. for 
ci es 


Mase 5, One has a mean miss distance of below 16m with 


} 


a Peete r ena 6 rad7see. The trend of the mean miss 
Meeeeaices 1s almost constant over certain range in values 
of es . The variation of the mean miss distance below 
mais ranece 1S a tenth of the magnitude comparing with the 
Beeemero) Of the mean miss distances as a function of 
Meeeamum acceleration. It is concluded that for scenario l, 
meemilean Miss distance in case 1 and case 2 is not effected 


Geeta oll] rate if the roll rate is above 3 rad/sec. In 


memo, Migher roll rate limit can reduce the mean miss 


Mewrance. For scenario 2, the mean miss distances of case 
lr Case 2 are jJike the trends for scenario 1. However, 


Miemmean Miss distance of case 5 has a minimum of 26m at 
ee mumeirad/sec, tien Ine€reases up to 36m then decreases 
fSeeeoilitly. {this trend may be due to the effect of the roll 
Meee ac, which is caused by the increased maximum roll 
meeenand by the underestimations of the state variables in 


Meeemo. lhe other characteristics, 1.¢., the mean miss 


Bemeence decreases with increased maximum roll rate limit 
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male tOr scenario |, case 2 shows less near miss distinc 


Meme case _, hut for scenario 2, the opposite trent she 
Meo latice with the Same reasons as those of tne Lirst 
Mia lsyvis. 

Bmem variation Of Lhe mean miss distances with the time 
Meemtant as a parameter is shown in Figures 5.5 and 5.6. 
Summer, Scenarios, one has ~a@ flat plateau for low values of 
Maeecime Constant and a sharp rise with increasing calues of 
Mme] CONStant, 1.¢€., for scenario 1, Figure 5.5 shows 


Mmepewecurves of the mean miss distances 


ih 


roi tinea nerease Wp 


45 


Memene time constant of 0.9 second, for scenario 2, th 


D 


Mees Up cO 0.5 second are plateau, then tor both scenarios 
Mmm@emcUrves rise sharply. For scenario 1, the mean miss 
Meeweamce Of Cass 2 is less than that or case 1 up to 0.3 
meeond time lag, but above 0.5 second lag, the results show 
Memeermesite trend. Aliso over this time lag, the curves of 
Memeo fOr bOth scenarios are not Sensitive to the variation 
Mmeearameter up to the time constant of 0.7 second for 
mueemario 1, 0.9 second for scenario 2. The general trend of 
Mm@emeurves, 1-€., as the time constant becomes longer, the 
ieee ciStance becomes larger, is due to slow system response. 
meemeetne figures showing the results, the time lag of less 
memo. 5 Second is necessary to obtain the mean miss 
feeweance Of the required order of magnitude. 


Pmeures 5./ and 5.8 show the variations of the mean miss 


MeeeaticeS aS a function of one sigma angle noise. For 
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Beeetario 1, as shown in Figure 5.°, the maxdffum mis 
Msecaitcs 15 ObtaINed For the Six mMeasureuen 


ir the two measurement 


CO 


Meees distance being the smailest £ 
Beppe DOve about the angie error of '.5 de 
Meedilce i5 NOt a strong function of the angie error. for 
Mem@@mmeie Or VaillleS, @-¢., less than 0.2 decree, the curves 
Meco be Qpproacnins each other. ter the scenario 2, the 
MeBE@ecistance 15 a Strong f£iimction of the angle error. For 
Mummers Oscetarios (Figure 5.) and $5.8), the two measilrement 
Mebeeeiias Similar results up to an error of 0.4 degree. Ail 
Meee S10h% a Strong rise in the miss distance with increas- 
Mrmeerre:. ihe ©TOllr measurement case seems pasticularly 
memsitive. 

[evomevyeariation of the miss distances as a function of 
Meme velocity are shown in Fistires 5.9 and 5.10. During 
Permeation the initial stand off distance was increased to 
meeoectic Cnoagement time the same. One important aspect of 
Mmeere 5.9 is that represents an opposite trend of the 
SememectiOns in the mean miss distance to the general trend 
Seeeene Variation in the mean miss distance. Under the 
Semeeee Velocity ot /00m/s, the mean miss distance is less 
Mapes ©tOr CaSe 3, from 15m to 12.5m for case 2, from 15./7m 
mame om £Or Case 1. However, the curves for scenario 2 
Show the usual trend, i.e., the mean miss distance of case l 


Memeene smallest and that of case 3 is the largest among the 


Seema ss Gistances Ot all three cases at a value of the 
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M@ereases, CxCept fOr One Curve of case 1. It may “ec sumna- 
Mmeeecd tie foregoing. Considerations as follows: As tne 
Memes Of the missile and target geometry components become 
Meee er, 1.¢6., relatively less efrect of neglecting the 
Beeeier Order terms in the first order fi 
Me@eamination of the measured vectors to noises, the mean 
eerecictance Can be enhanced with the estimator of the six- 
[Mermsourement vectors. AS the estimator algorithm becomes 
meme utile trend Of the mean miss distances is less sensi- 
[ee tO the target geometry. 

Memeo mast ieana)s. lJ °réeprcsent the Sifect on the mean 
MpeM@eCdceancc Gt che total ensagement time. For simulation 
Mm@ewtinictial stand off distance was increased to keep the 
Mme@eeant missile velocity. The general behavior of the miss 
Memeeeees Of all three cases slowly decrease as the engage- 
ment time increases. For case 1, the mean miss distance 
slowly decreases Or is almost constant at below 10m. In the 
Mmmmer two Cases, there is more variation in the miss distance, 
Mumeecehne ceneral trend is still downward with increasing time 


Memeo. Note that as time to go increases for a given T, 
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Mmeterc Wore Characteristic time interval availabie reflected 


ia |CL ois Pee tomcat ss apstance theretore should be less 
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femme cOtal enoaecment tame ancreases. This fact is veri- 


feed trOm both figures. 
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Merervacted its separate consideration. General trends of the 
Memes are Similar to the other scenarios but with higher 
Memeoecistance. ihis scenario is a strong test of the 
Meee. iin this scenario, the control laws are subject to 
MmemerroOr due to the limitation of the small anole assump- 


meemeand the higher maximum acceleration capability is 


tf) 
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ae 
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Guerre d aie tO the large variation of the PIC-mi 
meet ceulom. ihus there are some unexpected trends. 
Meet ee sels Shoxs the variation of the mean miss 


Seeeeances ior each case as the IiM@eteasces. ALchouch 


al 
THis 


M@enmissile veliccity is less than that of scenarios 1 anda 2, 


maemo ximum allowable acceleration is increased. As shown 


Memeerotre, the mean miss distance for case 1 is lln, 


meeemease | is 10m, for case 3 is 15.5m only at seen G Zoos 
Meme icure shows the miss distance variations in this 
Merar iO are more sensitive to the a Pei einetiie they 


Max 


Mememscenarios. the mean miss distance ot below 40m can be 
Mepeeemed at a = over 1199's, 1.€., more a... is 

max . max 
Mmeeered than in scenarios 1 and 2. This trend of the mean 


Geerecistances Can be explained with the larger variation of 
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angele 3%) Pavel rol eeriiand @euaAZ ions. 

Peevey 1) Fisure 3215, the variaticns of Epo, 
distances OCs ee mee oes tUnctron of 7p are imost 
Memstant above oe Pema sce, hac of case S, w0Wwever, 
@0ntinuousiv decreases up to vane ema SoG. «mei. 7s 
Meese CONStant. for scenario 3 Gone needs a maximum roil 


meee Of about 5 or 6 rad/sec to have an effective BTT 
Beeeele, which is higher than in scenarios 1 and 2. 

Meemre 5.15 shows the variations of the mean miss distances 
eee ct i0n Of time constant. Up to the time constant of 
Meese conc, the mean miss distance is less 720m. For all 
Meee casceS, COMparins with the results for scenarios 1 and 
Eyeene tigure shows the curves of the mean miss distances 
Meee SCNSIitive in this scenario. The target accelera- 
MuemercctOr lies in same direction as the missile velocity 
MemeorOwn in Figure 5.4. Because the components of the 

Meme ive velocity are major components in computing the 
Mmemterol commands, the more correct control commands may be 
Mememced due to the small uncertainty in Moy and Myo 
Semponents. This results in smooth variation of the mean 
Geeeewcistance curve. fhe trend of the miss distances as 
mmenone Sigma angle error increases, shows some different 
M@eeacteristics from the trend for scenarios 1 and 2 in 
Meee >. 16. The mean miss distances for all cases 


Menease CONtinuously as the noise input increases up to the 
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Siemsteia error Of anout (.+ degree. Above the ore sigma 


- 


then rises =u, 


Y} 


mreror Ol1 0.48 deyree , case 1 decrease 
meee COmvinves decredsins siouly, case 3 increases 
Gontinuously. The variations of the mean miss distances in 


SMeemario | are almost piateau and in scenario 2? are aimost 


Meistant Up tO the one Sigma angle error of about 0.4 degree 


Meter ises Siiarpliy as shown in Figures 5.;, and 5.8. This 
Mer onecan be dedticed from the fact that as the one sigma 


mere error becomes larece, the same magnitude of the two 
Meeee yvectOrs in diftrerent direction, i.e., relative 
eee GtTocetson and relative in Y-direction, having 
memeoame Order Of magnitude may be contaminated to the noise, 
mum Otm@er Mand, the other scenarios, only the components 
Memes -cirection being relatively large are contaminated 
Mmemcne noise. From this fact, it can be concluded that the 
Mec in a Sideé-engagcement case 1S more senSitive to the 
mmere error than the other scenario. 

Meee woyl/ Snows tne Variacions of the mean miss 
Meeetances as the missile velocity increases as a parameter 
Mmeeees HhOlding the total engagement time of 5 seconds and 
Merieine the relative position of the target). Up to the 
Meeeod te velocity of 700m/s the curves of the mean miss 


Mmeeemecse are PildteGau, attér that velocity the curves rise 


mmr y. Figure 5.8 shows the variations of the mean miss 





ewe eecsncos 2 fUNCctIOn Of the total engagement t rit ao. 


(eee eGetrosc™ OG) the missile and 


eS) 


noiding fic opale ae 
@areet and changing the relative pesition of the Werwet in 
Secnario). 

PieUEewo lo Soows ee mean Miss distances rise sharply 
Gver the engagement time of 6 seconds. This is somewhat 
Unexpected and mav be an artifact of the implementation 
eeecudre., ii) Other Scenarios the mean miss distance as a 
moet ron Of the missile velocity increases slightly as the 
Semele sveclOcity increases and that as a function of the 
Motal engagement time decreases slowly with increased 
Sm@eecement time. From above analysis one has a conclusion 
mie the Computed control commands should converge in the 
meet tame as the computed commands of the two- four- 
ieeeirement case shown in Figure &.32 and Figure 3.54 to 
obtain the low mean miss distance. However, it can be 
expected in this scenario that the computed control commands 
Geeenoc COnverse possibly due to the violation of the small 
Meee assumption in deriving the equations of the control 
Memmands {1l.e., in this scenario the Mo. component 1S 
Meeoee in Equation (5.7), that causes the large mean miss 
Meeotances}. However, in the case of relatively small magni- 


he 


Mm@ee Of the missile and target geometry components, t 
fyeeor Caused by the violation of smail angle assumption 1s 
meeeeeively small comparing with the errors caused by the 


contaminated O@mEemeoumoadastred vectors to noises, 
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MepeearvomsesOr Lae cCOmtrol aiputs, And the ¢stimati 
Pe eee CChOLSeU-GNenhiist-Order [ilter. This ract 
mericticameith the resules that the mean miss distance is 


Mmimeror case 1, l0m for Case 2, lim for case 5 at the 


ws 


eect ic velocity of 6)0m/S and the engagement time of 
feeenas, 25 Shown in Fieures 5.17 and 5.18. 

iitemre Mate figire 5.19 Shows the variations of the mean 
SPM Ge iecs 1Or Cace 5 Of all three scenarios as a func- 
Meeeor tie One Sigma angie rate error. The result repre- 
[ene ciat tune Variations of the mean miss distances due to 
Mivercased input values of one Sigma angle rate error are 
Mest GCONStant fOr scenarios 1 and 2. For scenario 3, the 
Mremiss CisStance Varies at Constant rate of inclination. 
[meres lOpe of the curve is relatively small comparing with 
the slopes of the mean miss Hee tance ena ves aS ene oOtner 
Memeemeters change. From this analysis, it mav be concluded 
Mupmemetiac GXposure of the measured vectors to the angle rate 


Meeorenolse 15 negievible in the case of present study. 
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VE. CONCLUSTONS 


Peete Ue nett eecclve Sidsed cuidance law and three 
Seeeetiaed halmay rilter equations have been derived from the 
Me@eimal COmtrol theory rOr a bank-to-turn missile with zero- 
mreeatd with 0.5 second time lag autopilots for nitch accel- 
Beeaeion and roll state. 

om. sero wonemothon are linearited around the present 
See itetion Of the missile with the assumption of small 
Mere Ol anetes.,  Wurins Computer Simulation, tne optimal 
memrrol laws and estimators were tested separately invoking 
miewcepardation theorem. In the simulation of the system 
Miteec ODtimal esStimatcrs were used, 1.¢€., estimator 
mee twO measturement vectors, with four-measurement vectors 
Mraewith six-measurement vectors. The system has been eval- 
Meee Or ciree Scenarios, tail chase engagement, head on 
smeecgement and side approach engagement. 

Meemimencediaiyses Of the Simulation results in Chapters 
Mmeereard 5, it Can be concluded as follows: The controi law 
Peewee ccssitiiy Simulated for a hypothetical bank-to-turn 
Peo wien pitch acceleration and roll rate autopilot 
femenin the imposed limits on three scenarios. Miss distances 
Meeeneccro-lage were negligible for all three scenarios 
@eeeow 0.5m). For the time lag case one had miss distance 


Memeo enor Scciarlo |, 0.8m for scenario 2 and 32.9m, with 
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te second jac and with ne! Pgs. bY itposing 
Peet fF emek ium 2 Of 25 On the scenario 3, the miss 
Pst amceewas rediticed tO below 5m miss distance. Cne would 
Meme ce a ich Ses requirement for a sid 

Mmmeteas Im the ©€ase Of Seenario 53. It is evident that the 
Mmesiie'’s tracking ability against the target is very 

meet cive LO the mMisStie ana target geometry, the missile's 
Meme vering Capabvility is aiso a very important factor. 
Seemal Control laws applied to the side approach engage- 
memo lace a Severe Strain on the system from the viewpoint 
Seeecibc Small angle approximation. Increasing the number of 
Bee urements did not result in increased performance in the 
Mere Ot Sidi ler miss distances for the missile except in 
Mmeetted Cases. Inis may be due to the increased ccn- 
Meeexity Of the system as the number of measurements was 
MMe (ecoeruUrthe;, in this study only first order extended 
Meets ere implemented. It is possible that better 
Mmemrormance cOuld be obtained at the expense of increased 


memoutational load, by using a second order extended Kalman 


maiter. 
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SPPENOTX 2) PROGRAM LTST IN 


ieee a Outme ol ISting of the computer program 
@yin the present study. Only one program for the six 


prem om. Case is provicea. Except for a few differences 





Meteo Ue eOUL Mo 81titer the procrams for the two other 


ee eee One te) procram included here. 
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